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Abstract: We study the Casimir force on a single surface immersed in an inhomoge- 
neous medium. Specifically we study the vacuum fluctuations of a scalar field with a 
spatially varying squared mass, + \/\{x — a) + V{x), where ^ is a smooth potential 
and A(x) is a unit-area function sharply peaked around x = 0. A(x — a) represents a 
semi-penetrable thin plate placed at x = a. In the limits {A(x— a) a), A 00} 

the scalar field obeys a Dirichlet boundary condition, = 0, at x = a. We formulate 
the problem in general and solve it in several approximations and specific cases. In 
all the cases we have studied we find that the Casimir force on the plate points in the 
direction opposite to the force on the quanta of 0: it pushes the plate toward higher 
potential, hence our use of the term buoyancy. We investigate Casimir buoyancy for 
weak, reflectionless, or smooth V{x), and for several explicitly solvable examples. In the 
semiclassical approximation, which seems to be quite useful and accurate, the Casimir 
buoyancy is a local function of V{a). We extend our analysis to the analogous problem 
in n-dimensions with n — 1 translational symmetries, where Casimir divergences become 
more severe. We also extend the analysis to non-zero temperatures. 
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1. Introduction 

The interaction of a quantum field with a material medium sometimes can be idealized 
by placing a boundary condition on the field at the interface with the medium. Then 
the effects of the medium can be interpreted as modifications of the zero point energy 
of the field due to the boundary condition on the surface. The classic example is the 
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force between two parallel, grounded conducting plates due to zero point fluctuations 
of the electromagnetic field first discovered by Casimir[0. Hence problems of this 
nature are known in general as Casimir problems. If the surface S is embedded in an 
inhomogeneous medium, then the quantum fluctuations of (p sense the inhomogeneities 
and give rise to a force per unit area on the surface even in the absence of a second 
surface. We will refer to this force as Casimir buoyancy. As "buoyancy" implies, the 
force is opposite to the force that acts on the quanta of the fluctuating field, at least 
in the cases we have been able to study. 

In this paper we examine the problem of Casimir buoyancy. We formulate and 
study the problem for the simplest possible cases. We consider a scalar quantum field, 
0, obeying the boundary conditions imposed by a sharply peaked background on an 
n — 1 dimensional hyperplane embedded in n-dimensional Euclidean space. We assume 
that the mass of the scalar field varies in the direction normal to the surface, i.e. we 
introduce an interaction ^V{x±)(j)'^ into the Lagrangian, 



The function A, normalized to / dxA{x) = 1, is assumed to be sharply peaked at 
x± = a. Usually a Dirac (5-function will do. We are particularly interested in the 
"Dirichlet limit" , where the coupling A goes to infinity and the field obeys the Dirichlet 
condition, = 0, on the hyperplane. It is easy to imagine problems to which such a 
formulation applies, where a quantum field is subject to forces on two different scales, 
forces at a high energy scale that can be idealized as a boundary condition, and forces of 
order the mass of (j) that can be regarded as a smoothly varying background.^ We believe 
that similar considerations apply to a gauge vector {e.g. electromagnetic) field in an 
inhomogeneous medium and to fermion fields. The generic problem of a quantum field 
constrained on a surface and modulated by other forces in the bulk also arises in brane 
world scenarios, where similar effects should also be expected. We have not considered 
buoyancy in a curved space-time where the effect could arise from an inhomogeneous 
curvature 0. 

As usual in Casimir physics, there is little intuition to guide us a priori. For 
example, there is no reason to expect the buoyancy force at a point x to depend only 
on V at the point x. In general we find that the buoyancy depends non-locally on V. 
However when the background field is smooth enough to admit a WKB approximation 
we find that the buoyancy reduces to a local function of the background field. Likewise 

^In the real Casimir effect the ideahzation of the influence of the metal or dielectric on the electro- 
magnetic field as a static background must be taken with some care. In particular some divergencies 
arising from this idealization could be absent if a dynamic description of the material is adopted B . 



c 





-2- 



we know of no argument to give the sign of the Casimir buoyancy. Should it be parallel 
to the force on the quanta of 0, or antiparallel? In all the examples we have been able 
to study we find that the force is opposite to the force on the quanta, hence a buoyancy. 

Casimir problems can suffer from divergences of two different kinds The first 
are the familiar divergences that afflict any quantum field theory. The Casimir energy 
for a fluctuating field in a time-independent background^ a{x) is the full one-loop 
effective energy, E[a], the sum over all one loop Feynman diagrams with arbitrary 
insertions of a, as shown in the first line of Fig. The low order Feynman diagrams 
diverge: the 1-point function diverges for n > 1, the two point function for n > 3, etc. 
As usual, these divergences are cancelled by counterterms that are polynomials in a, 
i^CT = ciO" -|- |c20"^ -|- . . . (some are shown in Fig. |l|). The Casimir energy becomes 
dependent on the renormalized parameters of the a-field dynamics. For example the 
need for the counterterm \c2(t'^ in three dimensions generates a dependence on the 
renormalized cr-mass. 

The second type of divergence is more interesting and more challenging. A re- 
alistic material medium cannot be idealized by a boundary condition at all energy 
scales. When the frequency of the fluctuating field is high compared to the natural 
scale of the interactions that characterize the material, its effects fade away. A back- 
ground that constrains all modes of a fluctuating field is unphysical, and can introduce 
divergences into physical observables that cannot be removed by standard renormal- 
ization methods]^. The origin of these divergences is quite clear from the Feynman 
diagrams. Each diagram involves integrals over the momenta carried by the external 
lines, <y{x) = J (ipe~*^^o"(p). If the background has a discontinuity in some derivative, 
say the k-th, then at large p, a{p) ~ p~'^ and the integrals diverge for sufficiently large 
number of dimensions. There are no renormalization counter terms available to cancel 
these divergences. They signify that the quantity under consideration, even a directly 
measurable one like the Casimir force, is sensitive to the high energy cutoffs, Q, above 
which the material no longer affects the field (say a{p) = for p > Q). 

In a given model there may be some quantities which admit a Casimir {i.e. bound- 
ary condition) description and others that do not. The boundary condition idealization 
shares with the effective field theory the notion of a separation of scales. The material 
structure is characterized by a high energy scale, the cutoff Q. Modes with energies 
below this scale obey a boundary condition, modes with energies at or above Q do not. 
If the boundary condition method is applicable, then physics at energies much lower 
than Q can be described by the boundary condition without reference to Q at all. For 
Casimir's original problem, parallel conducting plates, the plasma frequency c^piasma 

^For us a{x) — Xd{x — a) + V{x). 
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sets the high energy scale fl ~ ^piasma, and the plate separation, d, or rather hc/d, is 
the energy scale of physical interest. When cUpiasma ^ c/d the force between parallel 
plates is well described by the boundary condition calculation. It can be shown that 
the Casimir idealization works for Casimir forces between rigid bodies in vacuum and 
for local observables like the energy density outside the material (at distances greater 
than hc/Q), two examples of immense practical importance. Other observables are not 
so fortunate. It was shown in Ref. that the Casimir pressure on a sphere of radius 
R, for example, depends on the cutoff, so it is not possible to study the pressure, even 
when hc/R -C fl, without characterizing the material in detail. 

It is important to learn the circumstances under which there is an effective low 
energy description of buoyancy independent of the material cutoffs. We find that the 
answer to this question depends on the number of dimensions, n. In one dimension, 
where the surface is a point, there are no divergences of any kind and the Casimir 
buoyancy is independent of the cutoffs. Thus, for example, it goes to a finite limit 
as A oo. Higher dimensions can be studied relatively easily using the formalism 
developed in Refs. [|7|. For n < 2 the situation is the same as for n = 1. For 2 < n < 3, 
the Casimir buoyancy remains independent of the details of the surface, but the limit 
A — >■ cx) cannot be taken (it diverges like A""^). For n > 3 there is no separation of 
scales. The buoyancy depends on the details of the structure of the material. 

In the next section we describe the formulation of the Casimir buoyancy problem in 
one dimension. We recast it in terms of the Schrodinger equation Greens function with 
potential V{x), and express the buoyancy in terms of the bound and scattering states 
of V. We study both fixed A and the limit A ^ oo. In Section III we describe some 
important approximations and special cases: For deep and smooth V{x) we derive a 
WKB approximation; we study the case where V{x) + = (a turning point) which 
cannot be analyzed by means of WKB, proving buoyancy in this well. We study 

the Casimir buoyancy force in a thermal state finding the buoyancy is qualitatively not 
affected by a non-zero temperature. For a reflectionless potential we show that only the 
bound states matter. We study the buoyancy outside the range of V{x) and construct 
the first Born approximation as an example of non-local but simple result. In Section 
IV we go through some explicit, solvable examples. Then in Section V we generalize 
from = 1 to higher dimensions. 

2. Formulation of the problem 

The situation of interest is summarized by eq. ( |1.1|) . It might arise if a quantum field 
is coupled to one field characterized by a high mass scale that lives on the hyperplane 
x± = 0, and to another field V characterized by a lower energy scale. As usual we take 
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a(x) a(x) o(x) 
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Figure 1: The Casimir energy of a Euctuating field, (j), coupled to a time independent 
background field, a, via Cj = ^a{x)(j)'^ is proportional to the sum of all one loop diagrams. 
The sum must include the contributions of counter terms, polynomials in a, required to 
cancel the loop divergences. The structure of the counter terms depends on the number 
of dimensions, n. The second line in the figure shows the counter terms required in three 
dimensions where the 1- and 2- point functions are primitively divergent. 



V to be time independent and externally determined; we ignore the back reaction of 
(f) on V. The principal dynamical effects occur in the one non-trivial direction along 
which V varies, so we begin by studying the one-dimensional problem. 

2.1 General considerations 

In one dimension the Lagrangian for reads 

C = l^'- 1<P" - ^ K + AA(z - a) + V{z)) ct? (2.1) 
We are interested in the vacuum energy of the 0-field as a functional of A and V , 

a sum over the discrete spectrum, if any, and integral over the continuum, where the 
{uj"^} are the eigenfrequencies of the Schrodinger Hamiltonian, 

n = -(f/dx'^ + m^ + \A{x-a) + V{x), (2.3) 
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which will, in general, range over both discrete and continues values.^ So the problem 
can be recast in the form of a Schrodinger equation, 7i0 = E(j), with E = uj"^ = k"^ + m?. 
If we define Q{x,x',E) = {x'\j:^\x) to be the equivalent Schrodinger Greens function, 
then we can write the density of states, dN/dE, as 



dN 

~dE TT 



- Im y d'^xGix^x^E + ie) (2.4) 



where n is the number of spatial dimensions. The ie displaces the poles in Q into the 
lower half i?-plane, and insures that the propagator, ^(x, x', t) = J dEe~^^^Q{x, x', E + 
ie), is causal. Substituting for the density of states into eq. ( |2.2| ), the Casimir energy 
can be written 



£ = —lm [ dE^fE [ d''xg{x,x,E + ie), (2.5) 
2vr Jq J 

where we have assumed that the spectrum is positive definite*^ and we have introduced a 
cutoff Q above which the boundary must be characterized in greater detail. We already 
discussed in the Introduction the conditions (on n and a = AA + V) under which we 
can take Q —>■ oo without affecting the low energy physics. We will go over some these 
arguments again in Section ^. At the moment we assume all these conditions to be 
fulfilled so we can take the limit Q oo without placing further restrictions on a{x). 

2.2 Force on a sharp surface 

We are interested in the case where A(a; — a) is a Dirac 5-function, A(x) = 6{x), and 
A > since this repels the field (p from the "surface" x = a. If we restrict the analysis 
to 1 spatial dimension (the extension to n > 1 with translational symmetry in n — 1 
dimensions will be performed in Section 1) it is possible to express the Greens function 
Q in terms of the (simpler) Greens function, Qq{x' ,x, E) = \x) in the presence 

of V{x) alone. 



no = -^+m^ + V{x). (2.6) 
dx"^ 

First, we write the Lippmann-Schwinger equation^ for Q in terms of Qo and A(x), 



g{x,x',E) = go{x,x',E) - \ j dygo{x,y,E)A{y - a)g{y,x',E) 



(2.7) 



■^We set h = c = 1 until further notice. 

^Negative energy single particle states in an external field correspond to vacuum instabilities that 
we do not consider here. 

^Here we are considering a fully renormalized Hamiltonian Ti., a function of renormalized masses 
and couplings. 
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When A{x — a) —>■ 6{x — a) it is easy to solve for Q 
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G(r r' E)-G (r r' E) - ^ ^oi^^ E)go{a, x' , E) 
G{x,x,E)-Go{x,x,E) X ^^^^^^^^^^^^ (2-8) 

For the density of states, we require the integral over x of Q{x, x, E). Using the identity 

[ dxgo{x,a,E)go{a,x,E) = [ dx{x\ \a) {a\ \x) 

J J rto — /to — -c/ 

= ^So{a,a,E) (2.9) 

we obtain 

[ dxQ{x,x,E) = [ dxgo{x,x,E) - ———^ —^go{a,a,E) 

J J 1 + Xyo{a,a,E)dE 

dxGoix, x^E)-^ In (1 + A^?o(a, a, E)) (2.10) 
Substituting into eq. ( p. 51) , we find 

£ = -^lm[ dEVE^\n{l + Xgo{a,a,E + ie)) (2.11) 
2tc Jq dE 

where we have dropped a term of the form £q = 2^Im dEy/E J dxQolx, x,E + ie), 
which does not depend on a and therefore does not contribute to the Casimir force, 
—dE /da. Since we are only interested in the force, not the energy, we differentiate with 
respect to a, = —dEjda, 

Im/ dE^fE——hi{\^XQ^{a,a,E^%e)) (2.12) 
27r Jo oEoa 

where we have taken ^ oo since according to our previous discussions the limit exists 
and is finite. 

The analytic structure of the integrand of eq. (|2.12|) is important for our analysis 
(see Figure 0) because only the imaginary part contributes to JF, and the integrand is 
real for real E except at its singularities. The gives a branch cut running from 



® After this paper was completed Brian Winn, in a conversation with one of us (AS), pointed out 
that singularly perturbed Hamiltonians like ( |2.3| ) have been studied in chaotic billiards theory where 
they are called Seba billiards and generalizations of Eq. (2.S) can be found in the literature on this 
subject S. 
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-E = to oo, which we place along the negative real axis. For real, positive E there 
are two regions of interest. Above threshold for scattering, E > m? , Qq is complex, so 
the integrand has a cut with branch point at E = m? . For E real and below threshold 
^o(fl, CL-i E) is real, so the only singularities occur when the argument of the logarithm, 
1 + y^Qo{a,a,E + ze), vanishes, where the integrand has poles. When E is below the 
spectrum of TYq, Go{o,,a,E) is positive and 1 + XQo^a.a.E) cannot vanish. So there 
are no singularities in the domain Q < E < rn? unless Tio has bound states. Suppose, 
then, that Tio has bound states at Ei, E2, E3, ... Em- Because Qo ~ — E) near 

the j^^ bound state, is easy to see that 1 + XQQ{a,a, E) must vanish at a value of E 
between the j*^ and j + 1*"^ bound states of Hq. At each of these energies there is a 
contribution to the imaginary part from the ie in the argument of Qq. Therefore at 
least M — 1 poles contribute to JF. If the pole just above Em occurs at an energy below 

then there is one more. These contributions have a simple physical interpretation: 
They are the contributions to the Casimir force from the bound states in V{x) subject 
to the boundary condition Aip'{a) = Xilj{a). The analytic structure of the integrand of 
eq. ( p^.l2 ) in the complex ii^-plane is summarized in Fig. (Q). 

From these considerations it is clear that the problem is simplified if we rotate 
the integration contour to the negative imaginary axis. Making the obvious analogy 
to Feynman diagram methods, we refer to this as "Wick rotation" to the "Euclidean 
form" of the Casimir buoyancy. There is no contribution to the force from the semi- 
circle at large \E\, because forE^V we have Go{a, a, E) ~ Although this yields 
a contribution to £ (logarithmically divergent in the cutoff it is independent of a, 
and therefore does not affect the force. 

The result, 

1 r°° (9^ 
r=— dE^fE-—\n{l + \Go{a,a,-E)) , (2.13) 
27r Jq oEoa 

is particularly useful because the argument of the logarithm is positive definite. Eq. ( |2.13| ) 
can be integrated by parts without contributions at the limits because the surface term 
at cxD is independent of a and the surface term at £" = vanishes: 

^(a. A, V) = -^ [ dE^-^ In (1 + XQoia, a, -E)) . (2.14) 
47r Jq y/E oa 

Here we have restored some of the arguments on the function JF as a reminder of its 
important variation with position. A, and background field, V. It will be sometimes 
convenient to introduce the imaginary momentum k = y/E + w? (here E is the dummy 
variable in ( |2.14| )) in the Euclidean domain, 

T{a,\,V) = -— dKK ^ in(l + Ago(a,a,-V/t^-m^)) , (2.15) 

Jm \l K — TfX Oil \ / 
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Figure 2: Analytic structure of the integrand of eq. (2.12 ) in the complex E-plane. The left 

comes from scattering 



m 



hand cut comes from the v E, the right hand cut beginning at E 
states. The pole contributions are marked by solid circles. They lie just above the bound 
states ofV{x), marked with (g). 



This result can be rewritten usefully by introducing the Jost solutions to the 
Schrodinger equation, Tiip = —Eip (notice the minus sign in front of E, due to the 
Wick rotation), which are defined by their behavior as x — > ±00 0, 

lim ip^{K,x)e'^'^ = 1 

x— >oo 

lim V^-(K,x)e-''^ = 1 (2.16) 



The boundary conditions of eqs. ( |2.16| ) render the functions ifj analytic for Re k > 0. 



On the contour — 00 < ii^ < needed in eq. ( |2.15| ) the Jost solutions are real and fall 



exponentially in the limits given in eq. (|2.16| ). The Greens function can be written in 
terms of the Jost solutions, 

T{k) 



Qoix, y, -V -m?) = — — ?/'+(/€, x>)^ (ft;,a;<) (2.17) 

2k 

where x> (a;<) is the greater (lesser) of x or and T{k) is the transmission coefficient. 
Upon substituting for Qq in eq. ( |2.15D , we obtain 



^ a, \,V) = -— dKK — In 1 + — ^ 

27r Jo y/n^ - m? da \ 2k 



'ip^{K,a)'4j-{K,a)) (2.18) 



In the last equation we have restored the factors of h and c {E = — m? has units of 
1/length^ in these units). 
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2.3 A oo: the Dirichlet case 

When the strength, A, of the sharp background field becomes large compared to the 
eigenvalue, uj, the mode of the field with frequency u vanishes at x = a. If all the modes 
that contribute to the buoyancy have eigenvalues below A, then effectively the field itself 
obeys the boundary condition 0(a,t) = 0. As discussed in Ref. 0, this description is 
acceptable if the physical observable of interest remains becomes independent of A in 
the limit. In the case of buoyancy force in one dimension the limit A ^ cxd exists. 

Urn J^{a,X,V) = T^{a,V) = -— / -=—\ngo{a,a,-E) (2.19) 



which can be also written as 



^oo(a, V) 



1 

An 



dtp~^{n, a) 



Vk^ -m^ \^+{K,a) da 



+ 



dip (k, a) 



IK, a] 



da 



2.20) 



It is instructive to study the contribution to the Casimir force of the bound states 
in the A ^ oo limit. ^ Returning to the non-rotated formula eq. ( p. 121 ), taking A ^ oo 
and substituting the Jost solutions representation for the Greens function, we have 

1 Id 
^oo(a,V) = -— Im / dE—=—\n{ilj^{E,a)ilj-{E,a)). (2.21) 
ZTT Jq ^ E oa 

Remembering that Jost solutions are real for E < m^, it is clear that the only contri- 
butions in this range come from possible zeros in ip^{E,a), where the integrand picks 
up an imaginary part because the integration contour goes above the pole on the Re 
E axis. The locations of these zeros depends on a, so we label them as Ej±{a). Near 
the j^^ zero in ip^, 



ij^{E,a) = {E 



Ej±(a) + ie 



dip^{E,a) 



dE 



+ ... 



(2.22) 



where we have restored the ie to make the nature of the singularity clear. Substituting 
into eq. (|2.21|) and extracting the imaginary part, we obtain a very simple expression 



for the bound state contribution to J-'^, 



J'ooia,V)\ 



he 



bound states 



E 



dEj± (a) 



± a 



da 



d_ 

da 



he 



E 



EMa) 



(2.23) 



^The analysis for finite A is similar, but technically more complicated. In practice, following the 
discussion in the previous section, one has to find the poles of 1 + A Qo{a, a, E) say E(a), which are 
also eigenvalues of and proceed as before. 
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Here we have used the fact that the various partial derivatives are related by 



d^p = ^dE + ^da = (2.24) 
oE da 

along the contour ilj^{Ej±{a), a) = in the E, a plane. 

The simple form of eq. ( |2.23| ) has a simple origin: the { Ej±(a)} are the eigenen- 
ergies of in V{x) constrained to vanish at x = a, so the quantity in brackets is the 
bound state contribution to the Casimir energy for a scalar field in the potential V{x) 
forced to vanish at x = a. This is an equivalent formulation of the Casimir buoyancy 
in the Dirichlet (A — > oo) limit. From this interpretation it is clear that the number of 
terms in the sum in eq. ( p.23|) lies between — 1 and N, where N is the number of 



bound states in V{x). 

To complete this parametrization of the Casimir buoyancy we write the continuum 
contribution as an integral over the scattering momentum k, so another expression for 
the total buoyancy in the X oo limit, equivalent to eq. ( |2.20|) is 



d 



If the fc-integration were rotated to the positive imaginary axis, k — > in, two kinds of 
contributions would arise: a) from the cut from ^^■^."'^^^ and b) from poles in ip^ip~ 
corresponding to the same states counted in the sum over j. The pole contribution 
would exactly cancel the sum and the integral over k would yield eq. (|2.2(]| ) as it must. 



3. Approximations and Special Cases 

Casimir buoyancy is an unfamiliar phenomenon. Its properties are not readily appar- 
ent from the general expressions, eqs. ( p.20| ) and ( |2.25| ). In this section we study the 
buoyancy in special situations where it simplifies. First we look at smooth, deep poten- 
tials where the WKB approximation applies. Next we look at reflectionless potentials, 
where we show that only the bound states contribute to the buoyancy. Third, we look 
at the buoyancy in the domain beyond the range of V{x), and finally we study the first 
Born approximation, where the buoyancy is simple, but not local. 

3.1 WKB approximation 



If the potential is smooth and strong, a WKB expansion can be made for the Casimir 
buoyancy. We work with the Wick rotated form, eq. ( p.2(]| ). To apply WKB [1^ we 



introduce a fictitious into the Hamiltonian, H. "^^^ + W{x) (henceforth we 
will write W{x) for + V{x)) and write 

iP^iE, x) = exp (-iso(x, E) + si{x, E) + ... j . (3.1) 



and find 



Sq{x. 



E) = ± [ dy^-E + W{y) 



si(x, E) = -i ln{-E + W{x)) (3.2) 

so the domain E < corresponds to the WKB forbidden region (since W{x) > 0, which 
we have assumed from the outset) and gives real so,i,..., ip± and Qq. 

The criterion for the validity of the WKB expansion (with h set to unity), s[ <^ Sq, 
reduces to 

d 1 , ^ 

— ^=== « 1, 3.3 

dx ^yw{x) 

from which we conclude that W'{x) /W^^'^ should be small and in particular where W = 
+ V{x) should not vanish® [i.e. at turning points). Hence the WKB approximation 
applies to potentials which are "deep" in the sense that / dy^JW{y) ^ 1 and far from 
turning points where W{x) = 0. 

It is straightforward to construct the Greens function, in the WKB approximation, 

Go WKB(a, a, -E) = ^ = (3.4) 

2^yE + W{a) 

and, upon substituting into eq. ( |2.14|) , and performing the integral over E, we find 
ma)wKB = 1-^ I TT v/W^ + A ln(iy (a)) - ^/X^ - AW{a) In ^ " ~ 



d 

= -— fi(a, A)wKB- (3.5) 
oa 

is a local function of the potential, V{a). We call this function the quantum potential 
n{a,X). 

Equation ( p.5| ) simplifies considerably in the Dirichlet limit, 

^oc WKB(a) = ^ — A/m2 + V{a) = -^l^(a)wKB- (3.6) 



^If W ~ a;^" then W' /W"^^"^ ^ a; " ^ becomes arbitrarily large as x ^ 0. 
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A remarkably simple result for such a complex phenomenon. One general feature can 
be easily deduced from this form of the quantum potential Q: it decreases when the 
potential V{x) increases and hence the force felt by the plate is a buoyancy force. 
We will encounter this phenomenon throughout the rest of this paper in extensions 
(to regions where the simple WKB form is not valid, to non-zero temperature and to 
higher number of dimensions) and several examples. 

3.2 Points at which V{x) + = 

When the potential is smooth (with a definite length-scale b ~ ["jl^^]"^) but there 
exists an xq such that V{xq) = (or in general an interval where (V"(x) + m^)6^ <^ 
1) then WKB cannot be applied. The simple formula eq. ( p.6|) breaks down. For 
example in the case of symmetric potential + V{x) = x"^ the force must vanish at 
X = for symmetry, while eq. ( |3.6| ) predicts — > constant. 

To understand what really happens when the delta function reaches the point xq one 
needs a more clever guess than WKB. Moreover, adding more terms of the WKB series 
cannot help since the asymptotic nature of WKB means that adding terms improves 
the result for larger values of b'^{V + m^), so certainly not close to Xq. The formally 
correct procedure would be to find a differential equation 'similar' to the one we are 



studying but solvable and a smooth map from one to the other ||TT|. In this way we 
would obtain a uniform approximation near the turning point Xq. However such an 
analysis goes beyond the goal of this paper, so to clarify the situation we will assume 
that close to xq = we have (remember that W{x) > 0) 

W{x) = V{x) + - 00^"^ + O {x^) . (3.7) 

If we have 7^ W{0) = Wq <^ this constant could be reabsorbed by a shift of E 
and the range of integration in eq. (|2.24|) , we will consider this case at the end of this 



section. 

We then study the propagator in the neighborhood of xq = i.e. the propagator 
Qo of the equation 

-^"{x) + wV(/)(x) = E(j){x) (3.8) 

which is the usual harmonic oscillator problem. We now set = 1 and we will reintro- 
duce it only at the end of the calculation to have dimensionally correct results. 

From the two independent Jost solutions of this equation it is straightforward to 
write the the propagator for the harmonic oscillator as 



2-iE+i)/2 /I e;^ 
Qoia, a, E) = — F ) e^"^ H {a)HE^ (-a) (3.9) 



-,2 
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where H is the Hermite function, a generahzation of the Hermite polynomials related 



to the parabolic cylinder function [jT2[. Notice that the poles of the gamma function 
give the correct spectrum of the bound states. 

One can then calculate the buoyancy force for any A, by inserting eq. (|3.9|) into 
eq. ( p.l4| ) and performing the integral numerically. In the Dirichlet case A — oo however 
the results simplifies using eq. ( |2.19| ) to 

r„ = 4 r + £)(^2±SVi^ _ ^zHiSViW) _ 2a) . (3.10) 



In this case the expression for small a can be recovered by expanding the integrand 
in powers of a and integrating term by term (the expansion can be carried on to any 
order of a because the integrals over E converge): 




(3.11) 

which again exhibits the buoyancy phenomenon, ie. the Casimir force pushes the plate 
toward points at higher potential. Notice that the force vanishes as a — > xq = as it 
should from symmetry arguments.^ 

If V{xq) + w? = Wq and Wq ^ 1/6^ (6^ = l/u in this case) one can repeat the 
preceding derivation, change the lower limit of integration over E from to Wob"^ [Wq/uj 
in the case at hand) and write 

^ooh.o. ^^^'^VW^), (3.12) 

where we have defined 




m-l *7^(7^-2H. (3.13) 

and /(O) = 3.24... . Although we do not have an analytic expression for / it can be 
studied numerically with ease and to a great accuracy. In the opposite limit Wq/uj ^ 1 



^If we want to compare this result with an exact one we can choose n — A Poschl- Teller (see Section 
4.2) with TO = \/20 for which V{x) + m? ^ 20x^ near x = 0. In this case J^oxact /.^^approx — 0.94 as 



a <C 1. Such an agreement must be considered impressive, since the propagator is not a local function 
of the potential V and a local approximation to V does not guarantee at all a local approximation to 
G. Evidently, for sufhciently smooth potentials this turns out to be the case. 
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one can still shift the range of integration and then, by working out the asymptotic 
limit of the propagator eq. ( |3.9| ) as E/uj ^ 1, one finds the now familiar WKB result 
eq. ( p.6| ). We will apply eq. ( p.l2D to the case of Poschl- Teller potential in Section 

One can also study higher order zeros (or minima) of W , i.e. points where W{x) ~ 
x^" in a similar fashion. The Jost solutions of such a problem should be found, possibly 
in a series for small x, and by means of them the propagator can be written explicitly 
in the neighborhood of the minimum x = 0. Here we will not extend our analysis to 
those cases since no qualitatively new phenomenon arises. 

3.3 Temperature dependence in the WKB approximation 

One can also study the case where the field is in a non-zero temperature thermal state 
rather than a vacuum state. The Casimir energy is then the integral of the energy 



momentum tensor component Too evaluated on the thermal state. The result is ||T3 




\g^{a,a,E)). (3.14) 



One can Wick-rotate this expression as well, since the Matsubara poles are on the 
negative real E axis. After the rotation the Greens function becomes real (the i from 
\fE cancels with an i from coth) and one has an imaginary part from the cot function. 
Dropping a term oc AlogQ arising from the semicircle at = Vl one has 



= ^ dE^ Im| cot -''j]^ (1 + A6^o(a, a, -E)) . (3.15) 

In the limit e — one has Imcot(x — it) = 7r^^(5(x — nir). The final expression is 
then a sum over Matsubara frequencies. 



The study of this expression would require a paper on its own. Here we will limit 
our analysis of the WKB case, since a closed exact expression can be easily obtained in 
this case. To obtain it however, we have to go back to the original, non Wick-rotated 
WKB expression eq. ( |3.14| ) (we again use the notation W{x) = V{x) + > 0) 




Sx^-^lmj dEV E coth\ \ ^\og{l + — = | . (3.17) 

2^yE + ie — W{a) 
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The force in the Dirichlet hmit (the cutoff Q can again be sent to oo without 
problems) is 




J^^ ~ — ^Im f dE\fEcoi\i I (3^ 1 ^ ( ^— ) , (3.18) 

ATI da Jo y 2 j dE\E~W{a)+ie; ' ^ ' 

and by trading the E derivative with a W{a) derivative inside the integral, taking the 
imaginary part and the limit e ^ we find 



where the last expression defines again a local thermal quantum potential 

fi(a,/3) = -iv^coth(/5v^/2). (3.20) 



In the low temperature limit T ^ we have Vt —yW /A as it should (com- 
pare with eq. (|3.6| )). In the high temperature limit T — oo, /5 ^ one finds 
J^oo oc T"^dW / da. Contrary to what one could expect from the known 'classical limit' 



JF oc T/a of Casimir force between rigid bodies 11^ , p!3| , the force goes to zero when 
the temperature grows indefinitely. However here we are in a totally different limit, 
where the background potential is slowly varying and we are considering only the zero 
reflection term for Qq. The rigid bodies expansion of Q that gives rise to the well-known 
Casimir force actually is made up of non-local reflection contributions |TB . 



This expression is valid whenever WKB is valid, hence when the length-scale of 
the potential h is such that WlP' ^ 1. We have already discussed the buoyancy effect, 
i.e. the quantum potential has a maximum where V has a minimum. A non-zero 
temperature does not modify this prediction qualitatively as can be seen from Figure 



3.4 Reflect ionless potentials 

The continuum contribution to the Casimir buoyancy vanishes in the Dirichlet limit if 
a potential is reflectionless. For a reflectionless potential, ip^ik.x) ^^^yc*^^ as a; — > 
— oo, so ip~{k,x) = :j^^y^~^*{k,x). In this case ip^{k,a)ip~ {k,a) = \ilj~^{k,a)\'^/T*{k). 
Looking back to eq. ( p.25|) , we see that the imaginary part of the logarithm is indepen- 



dent of a so the force vanishes. Thus the Casimir buoyancy in a reflectionless potential 
is entirely determined by the bound states. 
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Figure 3: Turning on the temperature does not affect the Casimir buoyancy quahtatively. 
Here we plot Q for the potential W{x) = W+x^ and /? = 1, 0.6, 0.4 from up down respectively. 
The minimum ofWatx = is always a maximum of f^wKB ■ 



3.5 Beyond the range of the potential 

Here we consider a background potential which vanishes identically for |x| > b. The 
same results apply to a short range potential {V{x) oc e~^'^') in the limit x oo. For 
X < —b in the Euclidean domain, 



KX 



ip (k, x) = e 

i'^{^.x) = ^e--- + ^e-^ (3.21) 
T[k) T[k) 



The Greens function becomes 



go{a, a, - _ ^2) = — (l + R{K)e^^'') (3.22) 

Ik 

so the buoyancy can be written, 

H\a) = -^ r '^!l_ Un(l + ^{l + R{K)e^^^)] fora<-6. (3.23) 
27r Jm \J K^ -vn^aa \ 2k ^ ' ) 

For X > b a. similar expression holds with R ^ R, the reflection coefficient on the 
right. If the potential is symmetric, then R = R. If the mass, m, of the fluctuating 
field is non-zero, then JF falls like e~'^"^"' when a > b. The explicit form depends on the 
refiection coefficient in the Euclidean domain. In Section IV the explicit example of a 
5-function background is studied. 
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3.6 First Born approximation 

When the background potential is weak the Casimir buoyancy can be expanded in 
powers of V. The first term is simple when m = and A oo. A straightforward 
calculation gives the Greens function to 0{V), 

Goia, a,E) = ^ + ^ (^e^^^" £ dye-'''yV{y) + e-''"^ dye'^'yV{y)^ + 0{V^). 

(3.24) 

has no bound states at this order, so only the continuum contribution in eq. (p.25|) 
need be calculated. Straightforward evaluation leads to 

(where the integral is intended as the Cauchy principal value) which is a non-local 
functional of V known to mathematicians as the Hilbert transform H of [|^] at the 
point a. 



4. Examples 

In this section we report the Casimir buoyancy in three explicit sample background 
potentials and use them to study the domains of application of the approximations in 
the previous section. First we treat the potential i{i + l)/x'^ on the half-line a; > 0. 
Second we explore the family of Poschl- Teller potentials, —n{n + l)sech^x. Finally we 
study the 5-function, pS{x), which has been studied in other contexts [ll7[|. 

4.1 V{x) = £{e + l)/x^ 

The potential Vq/x'^ leads to a well defined problem on the half-line x > when Vq is 
positive. It is convenient to parameterize Vq by i{i + 1), so we can use results familiar 
from the study of three dimensional central potentials. Note, however, that i need not 
be integer. Since V{x) is positive definite there is no obstruction to taking m = 0, 
which we adopt to simplify the calculation. The formalism of Section II has to be 
modified slightly on account of the boundary at x = 0. In particular, the Jost solution 
ip~{x,E) has to be replaced by the solution, (f){x,E) oc je{kx), regular at x = 0. The 
other Jost solution is given by hf\kr). With this adaptation, and noting that the 
potential has no bound states, we can compute the Casimir buoyancy as an integral 
over scattering states, rotated to imaginary momentum, 
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where I = i+1/2, and and are modified Bessel functions. It is easy to verify that 
the ^ integral is convergent. In the Dirichlet limit the buoyancy is simply proportional 
to (for dimensional reasons it could not be different), 



lim JF(A, 

A— >oo 



he 



27ra2 



(4.2) 




Figure 4: Casimir buoyancy, —a? J- /he, for V{x) = l{i + l)/x^ in the Dirichlet hmit. The 
dashed curve is the exact result of eq. (pT^j. The solid curve is the WKB approximation, 
eq. p^j. 



The potential £{£+!)/ x"^ is neither weak nor refiectionless. It is however smooth and 
deep, so the WKB approximation should be accurate. J?-wkb(A, a) can be calculated 
easily from eq. ( |3.6[ ), and the limit A — > oo is particularly simple, 

hc/IjlTT) 

^WKB oo{£, a) = (4.3) 

The exact result and the WKB approximation are compared in the A — oo limit in 
Fig.|[ The difference is very small once ^ ^ 1. In fact WKB does even better: if we 
make the standard Langer replacement, + 1) ^ {i + |)0, [ill, then the WKB 

and exact results coincide within the widths of the lines in Fig. |. 

In Fig. 1^ we compare the exact result with the WKB approximation for finite A. 
Clearly WKB is an excellent approximation in this case. 
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Figure 5: (a) Ratio of Casimir buoyancy to the WKB approximation for V{x) = £{i + l)/x'^ . 
Up to down (red to blue) £ = 1, 3, 5; (b) Same as (a) with the "correction", £{£+1) {£ + lf- 

4.2 Poschl- Teller potentials 

The potentials V^^{x) = —n(n + l)sech^(x), n = 1,2, 3, are reflectionless and the as- 
sociated solutions to the Schrodinger equation can be expressed in terms of elementary 
functions. The Casimir buoyancy can be computed from the bound states alone. For 
simplicity we restrict ourselves to the X oo limit in this case. In practice, it is easier 
to compute the integrals in the Euclidean region using eq. ( p.2(]| ) since the scattering 

^''The apparent dimensional inconsistency in V^f"^ needs a word of explanation. If we begin with a 
dimensionally correct Hamiltonian, —cP/dx^ + w? — Vosech^(a;/6) and define a dimensionless unit of 
distance, z = x/h, then Vb6^ is the dimensionless potential which equals n{n + 1) in the Poschl- Teller 
problem. 
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wavef unctions, ilj^{E,a) are easy to construct but the roots of ilj^{Ej±{a),a) = are 
hard to find for n > 2. The necessary Jost solutions for imaginary momentum are given 
by 



^p^{K,x)= [ — tanhxj (-^ — 2tanha;j ... (-^ — ntanhx ) e*''^ 



(4.4) 



up to an inessential normalization. The necessary product of Jost solutions is given by 
V']^ (/t, a)V'r (f^iO,) = tanh^ a — 

'i/'^(/€, a)-ip2{K, a) = (1 — — 3fi;tanha — 3 tanh^ a)(l — + S/ttanha — 3 tanh^4)5) 

for example for n = 1 and n = 2. The required integrals can be performed explicitly 
for n = 1 and numerically for larger n. 

The case n = 1 is particularly simple and instructive: 

^PT/ X :,- tanha sech^a ^, 
Tf^im, a) = he — = (4.6) 

2 V — tanh^ a 

As usual the Casimir force is in the direction opposite to the gradient of V{x). The 
WKB approximation to the Poschl- Teller potentials is 

. ^ , tanha sech^a 

wkb("^, «) = ricn{n + 1) — , (4.7) 



m? — n{n + l)sech^a 

In order to understand how well the WKB result approximates the exact results we 
have plotted the force divided by \Jn{n + 1) for fixed /i = m/ ^Jrl{Jl^P[) = \f2. In 
these variables, the WKB result eq. ( [4.7| ) is independent of n and the the exact curves 
tend to the WKB curve as n is increased. This is no surprise since for large n the 
Poschl- Teller potential becomes more and more semiclassical. 

Even for large n, however, if yU ^ 1 we cannot use WKB since the minimum of the 
potential at x = is almost a turning point, i.e. iV{0) + m'^)lP' <^ 1 (remember we set 
6 = 1). In this case however we can resort to the harmonic oscillator approximation of 
Section p.2| which, for a — > and specializing to the Poschl- Teller potential, takes the 
form (the notation is the same as in Section |3.2|) 



47r \ ^n{n + 1) 



The results are plotted in Figure |^. As in the comparison with WKB the agreement is 
better the higher n, i.e. the more 'semiclassical' is the potential. 
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Figure 6: Casimir force for Poschl-Teller potential in the Dirichlet limit. Exact vs. WKB 
results. From down up (red, green and blue) we have the exact results for n = 1, 2, 3 and in 
on the top, in black, the WKB approximation, = m?/n{n + 1) = 2 for all the curves. 



4.3 5-function background 

The case of a localized background potential reaches the extreme limit when V{x) — > 
/3S{x). The Casimir force here is a non-local effect, as can be understood from the 
fact that the potential vanishes (almost) everywhere, and the overall sign of the force 
cannot be predicted a priori. The WKB approximation as developed in the previous 
section includes only the local semiclassical modification of the Green's function, which 
vanishes in this case. There is a semi-classical method for this case {e.g. a.s (3, X —>■ oo) 
which sums over classical paths that reflect from one surface to the other ||18|, |15| and 
accounts successfully for both the sign and the magnitude of the Casimir force. By 
applying the analysis of Section III.C with R{k) = R{k) = (3/{2ik — j3), we obtain 



J^{X,(3, m,a) = 



ndn d 



In 



2k \ 2k + 13 



(4.9) 



The special case (3 = \ has been studied in other contexts (see eq. (70) of Ref. [0), 
and our result agrees for that special case. For A 7^ /? however no new physics arises 
and we will not study that case here. More interesting is the limit A,/5 ^ 00 with 
m = in which case the integral can be performed yielding 



hcTT 
"24^' 



(4.10) 
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Figure 7: Casimir force for Poschl-Teller potential in the Dirichlet limit. Exact vs. the 
harmonic oscillator approximation, eq. ( |4.^ . On the vertical axis we plotted the ratio of the 
Casimir force and the harmonic oscillator approximation a ^ for n = 1,2,3,4. For all the 
points we chose such that V{0) + = 0, so x = is a turning point. 



the correct Casimir force for two impenetrable walls in one dimension at distance a 
from each other. 

Another interesting limit is A ^ oo, /9 — for m = which makes contact with 
the Born approximation of Section III.D, 



JF(A,/3,m, a)| 



m=0,A— >oo 



hc/3 
Ana 



(4.11) 



which agrees with the result obtained from eq. (|3.25| ). Here we have a non-local force 
which is still a buoyancy force. 



5. Beyond one dimension 

Field theories become more divergent in higher dimensions. Casimir effects are no 
exception, indeed they are more problematic because both the loop divergences and 
the sharp background divergences become worse. Methods for renormalizing the loop 
divergences, at least in dimensions where the interaction with the background field is, 
in fact, renormalizable, have been worked out in Refs. P, |l^, Divergences arising 
from the sharp background can be avoided by smoothing out the ^-function as described 
in Ref. H. 
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Once the buoyancy has been calculated in one dimension, the extension to higher 
dimensions can be constructed by the methods of Ref. 0. The core dynamics is the 
same in higher dimensions. The results can be summarized as follows: All the results for 
n = 1 generalize to any n < 2 without complication. For 2 < n < 3 the generalization 
succeeds except that the limit A ^ oo cannot be taken. In that range of n the Casimir 
buoyancy depends on the cutoff, A, on the strength of the boundary interaction. Thus 
there is no separation of scales, no "effective" low energy description of Casimir buoy- 
ancy. As n — i> 3 two further complications arise: first, a new counterterm, Cqt = C2cr^, 
must be introduced to renormalize the two point (in the background field, a) func- 
tion. This means that the buoyancy depends on a renormalized coupling, actually the 
mass of the a, that has to be specified. Second, as suggested by the appearance of 
an interaction proportional to cr^(x) induced by renormalization, the ^-function gives a 
divergent buoyancy. Instead it is necessary to smooth it out, replacing it for example, 
by a Gaussian, as described in detail in Ref. 0. The buoyancy force then depends 
explicitly on the structure of the surface as well as the shape of the background. In the 
case of a real material one should make oneself sure of giving a proper description of the 
reaction of the material, using for example a plasma model for the metal with plasma 
frequency cUpiasma- In that case we expect that, if the A — oo and the cUpiasma — ^ oo 
limits exist then they coincide. 

The problem of interest is an n — 1 dimensional hyperplane immersed in a medium 
which is modulated in the x_L-direction. As shown in Ref. 0, the Casimir energy (per 
unit 'area' of the n — 1 transverse directions) for such an interface can be written as, 

(5.1) 

where \fE = y/k"^ + w? is the contribution to the energy from momentum in the 
x± direction. dN/dE is the density of states for the one dimensional problem and 
dN^^^ /dE is the contribution to the density of states to M^^ order in the background 
field, (t(x) = AA(a; — a) + V{x)}^ £ft){M) is the contribution of the Feynman dia- 
grams through M^^ order in the background field plus the contributions of counterterms 
necessary for renormalization. 

The number of Born subtractions is determined by the degree of divergence of the 
field theory in a with a ct0^ coupling. For n < 3 only the tadpole diagram diverges, so 

^^The "Levinson" subtraction discussed at length in Ref. 0, is unnecessary here since the first Born 
subtraction does not contribute. To carry through the analysis for n = 2 a Levinson subtraction would 
be necessary. 
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only the first Born approximation must be subtracted. For the moment we set = 1 
and discuss the extension to n = 3 at the end of this section. 

The p-integration in eq. { f).l\j is to be understood in the sense of dimensional 
regularization, and can be performed/^ 

The first Born approximation and the contribution of the tadpole diagram plus coun- 
terterm are linear in A and linear in V . Therefore the V dependent term they generate 
is independent of a, so they do not contribute to the buoyancy and can be dropped, 

ir(-rV2) r,r,r,n/2d'N 

This remarkably simple form is valid for n < 3, and, of course, agrees with Section II 
for n = 1. 

As in one dimension, we take o{x) = X6{x — a) + V, construct the Greens function, 
rotate the integration to the negative E axis and integrate by parts, 

1 rn _ n) sin BE f) 



This should be compared with eq. ( 2.15 ), which is the specialization to n = 1. 

Note that the poles in r(l — n/2) are cancelled by the factor of sinn7r/2, which 
arose from the imaginary part of (— i?)"/^. Thus the only singularities in eq. ( ^.4|) arise 
from divergences at the upper limit of the E integration. These get more serious as 
n increases because of the prefactor of E'"/^"^. It is easy to see that the E'-integral in 
eq. ( |5.4| ) converges for n < 3: The high (negative) energy behavior of G{a, a, —E) is 
dictated by the WKB approximation, eq. (K 



G(a,a,-E)^ — + (E~^) (5.5) 

(here again W = + V) and this behavior leads to a convergent integral for < 3. 
However, if we attempt to take A — *■ oo before doing the integral, we lose convergence 
already at = 2. 



^^Formally, the integral is defined and performed for values of Re n small enough that it converges. 
It is then analytically continued to real, positive n, with careful treatment of singularities encountered 
along the way. 
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Beyond eq. ( ^.41) the buoyancy depends on the specific form of V{x). To get a feel 
for its variation, we evaluate the WKB approximation, 

^KB = ^A^W-^'Zr f 1 - - (nl) f ( ^4^] , (5.6) 



2(47r)t+i da \ 2) \ 2) ' \2^W{ 

where 

Jo y + 1 + 1 + X 

This function can be expressed by means of hypergeometric functions as 



m = --({l~x')^ -l)csc(-)+—Ti---)T(-)F(l 



X ' 'J ^"^ \ 2 J ' \2 2)' \2) ^ V" 2 2' 2' 

(5.8) 

however does not carry much more information than the original integral, in particular 
because the interesting limit A — oo translates in a; — >■ cx) where the hypergeometric 
function has both a real and an imaginary part. The imaginary part is cancelled by the 
first term and the real part is the one we are interested in but this is difficult to isolate. 
However, before carrying on an asymptotic analysis of this integral it is straightforward 
to evaluate the n = 2 result 



where the quantum potential is 



a- ^ . (^v - ^) iog(i + ^) - w ,og,2fl) I . (5.10) 

Neither the quantum potential nor the force have a finite limit when A ^ oo in 2 
dimensions. They diverge logarithmically in A. The case n = 1 result coincides with 
the result of Section II. 

Let us now return to the original integral representation to get an asymptotic 
expansion for large x {i.e. large A). Writing it as 



noo POO -| 

/ dze-'''2 / rf?/e-^^-(i/2 - l)"/2-i, (5.11) 
Jo Ji y 

performing the integral over y, expanding in series for small z and integrating term by 

term in z one gets^^ 

1 1 fmi\ 1 

/ [x) ~ 2-K csc(n7r) H — vr esc ( — j + (2 — n)TX csc(n7r) — 



^•^Alternatively one could use the Mellin representation of the hypergeometric function | [T^ and move 
the contour of integration in order to pick the desired poles. This procedures is much more involved 
than the one described here and yields the same results. 
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This expansion agrees very well with the exact expression eq. ( [5. 81 ) already at x ~ 2, 
as can be seen in Figure §, except in the limit n — 1. The expansion ( ^.12 ) has a pole 
in n = 1 that goes like x~^. This is not present in the original function / and indeed 
adding more terms pushes the pole to higher and higher powers of x~^. However in this 
case we know the exact result either from Section II (equation ( p.5|) ) or by taking the 
limit n ^ 1 of eq. ( ^.8| ) before expanding for A ^ oo (the two limits do not commute). 

One thing to notice is that in passing from n < 2 to n > 2 the first term in 
eq. ( 5.12|) becomes dominant over the second one. At exactly n = 2 they almost cancel 
each other leaving behind a term log(a:)/x which reproduces the leading term in eq. 
eq. ( |5.9| ) (also the higher order terms agree with the expansion of eq. ( ^.9] ) for A oo). 
Another thing is that from this asymptotic expansion one can see that in the Dirichlet 
limit one has a finite force if and only if n < 2. For n < 2 the second term in eq. (|5.12|) 
is dominant over the first one and hence the force is independent of A 

^«^^Re^^,y»?r(l-^). (5.13) 

°° (ia4(47r)2 V 2/ 

Having neglected the first term in eq. ( ^.12|) , however, this form cannot be used when 
n ^ 2. Incidentally however notice that the limit n — 1 now can be taken safely (and 
it gives the usual result eq. (|3.6| )) since we have not included the term oc csc(n7r)/x^~" 
which generated pole in n = 1. This pole was a fiction of our procedure of taking the 
limit A — s> oo before n ^ 1. It was signifying that for n = 1 the next term in the large 
X expansion of / decreases slower than l/x^, indeed it is oc log(x)/x^. As we said, 
including more and more terms in the asymptotic expansion pushes this pole further 
and further away. 

As anticipated, for n > 2 the buoyancy force diverges when A ^ oo. Keeping the 
first term in eq. (|5.12| ) one obtains 

^r-^?^A-rfl-?)secf4). (5.14, 



4'^7r2 da \ 2 J \ 2 

Notice that for n — > 3 we have a structure JF oc ■:^^XdV/da which comes from a 
term cx log A J d'^x{V (x) + X6 {x — a)y in the effective action (here A is the QFT cutoff). 
The pole l/(n — 3) from sec(n7r/2) is the dimensional regularization way of seeing the 
logarithmic divergence log A in the two- legs graphs for 3 + 1 dimensions. 

The result eq. (|5.6|) (and its asymptotic expansions eq. ( |5.13|) and eq. (|5.14| )) em- 
bodies all the behaviors we expect from the general analysis. It is a continuous function 
of n and finite for n < 3 and the Dirichlet limit, A — > oo, can only be taken for n < 2. 
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Figure 8: Comparison of f in eq. (5.8) with its asymptotic expansion eq. (5.12) for various 
n from up to down n = 1.4, 1.8, 2, 2.2, 2.6. 



To go all the way to n = 3, the plane in three dimensional space, it would be 
necessary to invoke the full apparatus of the field theoretic approach to Casimir effects: 
the surface must be smoothed out and the second Born approximation to the density 
of states must be subtracted from the ^'-integration and added back together with the 
counterterm ClS db Feynman diagram. 



6. Is Casimir Buoyancy universal? 

In the previous examples we have found Casimir buoyancy is ubiquitous. However 
all the previous examples share the fact that they can be very well described by the 
zero-order reflection term in the propagator. The WKB approximation was always an 
excellent approximation. There are obvious examples in which this is not true: Take 
for example the background V as made of a tall wall (a tall square wall will do the job) 
plus a shallow potential, decreasing toward the wall. If the shallow potential would not 
be present the delta function would be attracted toward the wall so for continuity this 
must be true for sufficiently small potentials even if the total potential is decreasing 
toward the wall. This is true because in the expression for the propagator Qq we cannot 
neglect the influence of the tall wall (that is why we called it 'tall') and this enters only 
through a non perturbative term in the form of the contribution from a closed path 
that goes from a to the wall and bounces back to a. More quantitatively by considering 
the background as made up of a wall and a smooth W, one can expand the propagator 
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as |2g 



-«7r/i/2 



Qoia, a, E) = , ' + " e^^^^"'"'^) (6.1) 

2^E - W{a) 2^E-W{a) 

where S'i(a, a, E) is the action of the closed path that goes from a to a bouncing on the 
wall and /i is the Maslov index of this orbit (/x = 1 is the wall is 'smooth' on frequencies 
\/E and /i = 2 for otherwise 'hard' walls). Inserting into eq. ( p.l4[ ) one finds a general 
formula including both local (e.g. W and its derivatives) and non-local (depending on 
integrals of W or independent of W) contributions. The two effects interplay and they 
are difficult to separate but it is easy to extract the two limits: 1) When — > one 
recovers the Casimir attraction between impenetrable plates (compare with eq. ( |4.9|) 
with (3 — >• oo) 2) When the wall is far away {Si 3> 1) the second term is exponentially 
smaller than the first and we find buoyancy from the first term only. 

In general then a more correct statement would be that Casimir buoyancy will occur 
in the cases where non-local effects, like closed orbits contributions to the propagator, 
are negligible. We believe we have presented enough examples here to convince the 
reader that this is not a too stringent request. 
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